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' ABSTRACT 



A recent laboratory experiment (Blaauwgeers et al. 2003) suggests that a Kelvin- 
Helmholtz (KH) instability at the interface between two superfluids, one rotating 
^Cgj ' and anisotropic, the other stationary and isotropic, may trigger sudden spin-up of 

I/" - ) , the stationary superfluid. This result suggests that a KH instability at the crust-core 

IO ' ( 1 S'o- 3 P2-superfluid) boundary of a neutron star may provide a trigger mechanism for 

pulsar glitches. We calculate the dispersion relation of the KH instability involving 
two different superfluids including the normal fluid components and their effects on 
stability, particularly entropy transport. We show that an entropy difference between 
the core and crust superfluids reduces the threshold differential shear velocity and 
threshold crust-core density ratio. We evaluate the wavelength of maximum growth of 
O ' the instability for neutron star parameters and find the resultant circulation transfer 

+3 . to be within the range observed in pulsar glitches. 

c/3 . 
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interiors - stars: neutron 



> 
ON 



in 
o 



1 INTRODUCTION 

It has been proposed that the neutrons making up the inner crust and core of a neutron star display superfiuid-like behaviour 
(Tamagaki 1970; Yakovlev, Levenfish & Shibanov 1999). They are compressed to such a high density that they act like 
composite bosons, forming Bardeen-Cooper-Schrieffer (BCS) pairs which condense into the ground state below a critical 
temperature T c and behave like terrestrial superfluids (e.g. 4 He). Calculations of the pair interaction and energy gap suggest 
that different types of Cooper pairs populate the crust and the core of a neutron star, 3 P2 (triplet) pairs in the core and 1 So 
(singlet) pairs in the crust (Hoffberg et al. 1970; Tamagaki 1970; Baym & Pethick 1975; Brand & Pleiner 1981; Yakovlev, 
Levenfish & Shibanov 1999; Khodel, Khodel & Clark 2001). Superfluidity enhances neutron star cooling, reduces the heat 
capacity of the star, and suppresses neutrino emission (Yakovlev, Levenfish & Shibanov 1999; Tsuruta et al. 2001). In addition, 
pulsar glitches (discontinuous spin-up events) are attributed to the motion of superfluid vortex lines; these lines tend to be 
pinned to the nuclei in the stellar crust and sudden, large-scale creep of these lines from one pinning site to another may be 
responsible for glitches (Link, Epstein & Baym 1993) 

Macroscopic phase coherence endows a superfluid with frictionless flow and irrotationality (Putterman 1974; Tilley & 
Tilley 1994; Khalatnikov 2000; Vollhardt & Wolfle 2002). Irrotationality leads to formation of discrete quantized vortices 
if the fluid rotates. This behaviour is well known and observed in terrestrial samples of superfluid 4 He (Putterman 1974; 
Tilley & Tilley 1994). Landau (1941) put forth what is known as the two-fluid model (where 4 He is modelled as two separate 
components, the inviscid superfluid and the viscous normal fluid) to explain the various phenomena observed in superfluid 
flows such as second sound (an entropy-temperature wave, where the two components move in opposite directions, induced by 
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a temperature gradient) and the fountain effect (Putterman 1974; Khalatnikov 2000). This hydrodynamic model is successful 
and has been polished to a high degree of sophistication, taking into account more and more effects, such as the motion of 
quantized vortices and the drag forces they exert (Hall & Vinen 1956; Barenghi & Jones 1988). Superfluidity in 3 He, on the 
other hand, was first observed more recently (Leggett 1975; Wheatley 1975) (see also Lee 1997 for a first-person account of 
its discovery) and the search for a complete hydrodynamic theory continues; for a comprehensive and up-to-date account, see 
Vollhardt & Woffle (2002). 

3 He superfluid is a condensate of 3 He Cooper pairs. The generally accepted theory of 3 He superfluidity postulates p-wave 
pairing with spin angular momentum h and orbital angular momentum h. Experiments show that there are two distinct phases 
of superfluid 3 He, known as He-A and He-B. He-B occupies the 3 Po state; it is isotropic and behaves like 4 He superfluid. 
He-A occupies the 3 P 2 state; the axis of the orbital angular momentum defines a preferred direction (generally denoted as 
1), breaking the symmetry of the ground state. 1 This anisotropy manifests itself as a 'texture', defined to be the vector field 
tangent to the local direction of the pair angular momentum (Leggett 1975; Lee 1997; Vollhardt & Wolfle 2002). The texture in 
the lowest energy state (undisturbed He-A) is uniform, but heat and fluid flows, applied electromagnetic fields, and container 
walls deform it. The texture, in turn, inhibits flows and heat transport perpendicular to itself. Hence, the presence of textures 
modifies the characteristics of any flows induced upon the fluid. 

A recent laboratory experiment with superposed He-A and He-B in differential rotation (at an interface stabilised by 
an applied magnetic field) indicates that the Kelvin-Helmholtz (KH) instability at the interface plays a significant role in 
triggering transfer of circulation from He-A to the initially irrotational He-B, spinning up the He-B sample; this is detected as 
discrete jumps in nuclear magnetic resonance absorption signals (Blaauwgeers et al. 2003; Volovik 2002). This exciting result 
suggests that a similar process may occur at the crust-core superfluid interface in a neutron star, with important implications 
for pulsar glitches. 

In Section 2 we present an analysis of the KH instability involving two superposed superfluids, taking as our starting 
point the two-fluid model, modified for 3 He. Our analysis extends earlier work in two ways: we treat the flow of entropy in 
the system self-consistently (cf. Volovik 2002) and consider two different superposed fluids (each with its own superfluid and 
normal fluid components) in relative motion, not the superfluid and normal fluid components of one fluid (Korshunov 2002). 
In section 3, we apply our KH dispersion relation to the experiment conducted by Blaauwgeers et al. (2003) and to an idealised 
neutron star model. This work complements an analysis by Andersson, Comer & Prix (2003) of the two-stream instability, 
a different process involving two interpenetrating, isotropic superfluids, which is also a possible triggering mechanism for 
glitches. 



2 KELVIN-HELMHOLTZ INSTABILITY 
2.1 Two-fluid model of 3 He 

We restrict our attention to the hydrodynamic regime, where disturbances of the equilibrium are much slower than the typical 
particle collision time. The nondissipative hydrodynamic equations describing an anisotropic superfluid are (Putterman 1974; 
Landau & Lifshitz 1987; Vollhardt & Wolfle 2002) 

( dv n ,i \ p n dp dT psP-a , „ f d 2 z' n d 2 z'„\ 

Pn + v„ • Vv n>i J =--—-sp s —- — (v„ • VK S , + Pn g, + a nti + J 6{z), (1) 

( dv s ,i \ p s h t dl k dl m p s dp dT p s p n ( d 2 z' 3 d 2 z' s \ . . . 



\ 8t ' ' s y " 3 y 2m^ Km " J dt dxi pdxi' "dxi' p y n3 -■*\ k ai ; a 1 dy 

§1 + Vn • v) ps = 0, (3) 



+ v n • V ) k = -\4.-^£ijklj(v„ s ■ V)Z fc , (4) 



V • v„ = 0, (5) 

V • v s = 0, (6) 



g- t + v n • V ) p n = 0, (7) 

l+v.-V)/) =0. (S) 

1 The experimental identification of the He-A and He-B as anisotropic and isotropic Cooper pair superfluids remains ambiguous (Wheat- 
ley 1975; Tang et al. 1991). 
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(10) 



The interface is defined to be at z — 0, the zeroth order flows are in the :r-direction, and the velocities in (9) and (10) 
are evaluated at the unperturbed interface. In (1)-(10), p s , n are the super- and normal fluid densities respectively, g is the 
gravitational acceleration (directed towards negative z-axis), er„, s is the surface tension force per unit length for the normal and 
superfluid interfaces respectively (directed towards the positive z-axis), z' niS is the perturbed z-coordinate of the normal and 
superfiuid interfaces, p, T, s, and m denote the pressure, temperature, specific entropy, and mass of an individual constituent 
particle, 1 is the unit vector tangent to the texture (in other words, along the direction of the pair angular momentum), v s ,„ 
are the superfluid and normal fluid velocities, and v„ 3 = v s — v n is the counterflow velocity. Note that there is no unique 
equation of state relating p s , p n , s,T, and p; these quantities evolve independently through the heat transfer (energy) equation 
(3). An isotropic superfluid is a special case, in which one sets 1 = in (1)-(10). 

Equations (1) and (2), the equations of motion governing the normal fluid and the superfluid, are modified with respect 
to a classical fluid by adding the terms sp s dT/dxi and (p s pn/ p)(v„ s • V)«„ Sj i. These terms come from the gradient of the 
superfiuid chemical potential p. In a classical fluid, one simply has V/i = Vp, whereas in a superfiuid, the entropy gradient 
can drive flows separately from the pressure gradient, because entropy is carried only by the normal fluid component. The 
sp a dT/dxi term leads to effects like second sound, a temperature-entropy wave which differs from 'first sound', the usual 
pressure-density wave (Landau 1941). In addition, the first term on the right-hand side of (2) describes the force exerted 
by the texture in an anisotropic superfluid when the texture is bent; the ground state of the superfluid far from any walls 
corresponds to uniform 1 (see Section 2.2), so the texture 'stiffens' the system and opposes flows which might bend it (Vollhardt 
& Wolfle 2002). Equation (3) is our (non-dissipative) equation of entropy transport [as for the electron-and-muon fluid in 
Mendell (1991)]. Equation (4) describes the evolution of the texture. Bending the texture (i.e. changing the direction of I) 
shifts the energy gap and the momentum distribution of normal fluid excitations, which then reacts back on the texture 
itself. The term on the right-hand side of (4) derives from the dependence of the free energy density (which is a function 
of dh/dxj) on the superfluid velocity (Vollhardt & Wolfle 2002). The remaining equations are standard: (5) and (6) enforce 
incompressibility, (7) and (8) enforce mass continuity, and (9) and (10) ensure that no fluid particles cross the interface, i.e. 
the fluids are immiscible (Chandrasekhar 1961). 

For simplicity, we make some approximations. We set all coefficients of conductivity on the right-hand side of (3) to 
zero. We also neglect the anisotropy in density for flows along and transverse to the texture, and we omit terms involving 
changes in free energy as the local texture changes, i.e. we postulate that the texture is in thermodynamic equilibrium to 
the zeroth order. In reality, (2) features a heat current term proportional to VT, with the thermal conductivities for flows 
parallel and perpendicular to 1 scaling as k« ~ (T c /T)kn(T c ) and n± ~ (T/T c )kn(T c ) respectively, where kjv(T c ) is the 
(isotropic) conductivity at the critical temperature. Similarly, although the physical density of the fluid at any point is 
uniquely defined, the anisotropic texture modifies the mass current to become p±v perpendicular to 1 and p\\v parallel to 1. 
This counterintuitive property, which arises because the two-fluid description is an effective theory where it is not possible 
to label individual 3 He atoms as s or n, is more pronounced at lower temperatures, with p n y = 7r 2 (m*/m)(fcsT/Ao) 2 p and 
pn± = (77r 4 /15)(m*/m)(fc_gT/Ao) 4 p as T — > 0, for example. Here, m* is the effective mass of each constituent particle and 
Ao is the maximum gap parameter, the energy needed to excite one pair out of the ground state (Vollhardt & Wolfle 2002). 



2.2 Unperturbed interface: uniform texture 

We investigate the situation where an isotropic superfluid is superposed on an anisotropic superfluid and moves relative to 
it with a constant and uniform velocity. We follow the derivation of the dispersion relation of the classical KH instability 
described by Chandrasekhar (1961), assuming continuous stratification of the fluid variables in the top and bottom layers 
with a discontinuous jump at the interface. The normal and superfiuid components interpenetrate; their surfaces are allowed 
to oscillate with different amplitudes, but must be in phase with each other. 

It is clear that a situation where all the quantities are uniform and constant in time satisfies (1)-(10). We take this 
equilibrium solution as our starting point. A uniform and constant texture, along with uniform and constant normal fluid 
velocity, satisfies (4) identically. A nonuniform but constant texture is also possible. In fact, it can be shown that the most 
energetically favourable texture for a superfluid in a spherical container is one where the texture lines radiate from a singularity 
on the wall, called a 'bouquet', 'fountain', or 'boojum' (Vollhardt & Wolfle 2002). In this paper we assume a uniform and 
constant texture without any singularities, as well as uniform and constant fluid velocities. 
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2.3 Perturbed interface: surface Kelvin-Helmholtz wave 

We linearise (1)-(10), and take all perturbed (first order) quantities to be proportional to exp(ik x x + ik v y — iwt), assuming 
that the wavelength of the perturbation is much smaller than the length-scale over which the texture changes. We take all 
zeroth order quantities (except temperature, which is the same for both fluids) to be constant and uniform except in the 
z direction, allowing for vertical stratification. The zeroth-order shear velocity is chosen in the x direction. From (2), one 
can see that, given a spatially uniform zeroth-order texture that is also constant in time, the textural contribution to the 
momentum flux is of quadratic order and thus negligible. One can now ignore the textural effects on the fluid flows, although 
the first-order texture still evolves according to (4). Writing D = d/dz, we then obtain 

~~ P" fAU \J „0 0/-, x PsPn 



ip n {-u + k x U n )u n + p n {DU n )w n = Pn J ik x ) P '-s°p°{ik x )T- Z n [{U s -U n )ik x (u s -u n ) + (w s -w n )D{U s -U n )UU) 

Pn + Ps Pa + Pn 

ip n (-u) + k x U„)v n = Pn (iky)p' - s p° a (iky)r PaP " (U s - U„)ik x (v a - v„), (12) 

Pn + Ps Pn + Ps 



S(Z), (13) 



ip n (-u) + k x U n )w n = Pn Dp - s°p1Dt - gp' n ( f s f" n (U s - U n )ik x (w s - w„) - k 2 a n z' n 

Pn + Ps Pn + Ps 

ip° s (~LU + k x Us)u s +p° s (DUs)w s = ~f S (ik x ) P ' + S°p°(ik x )T + f S f n [{Us ~ Un)ik x (u s - U n ) + {w s - W n )D(U s ~ U n )],(U) 

Pn ~r Ps Ps "T Pn 

ipi{~io + k x U s )v s = Q _f 5 (iky)P + s°p ( s(ik y )r + E*£lL- {jj a - U n )ik x (v s - v n ), (15) 

Pn * Ps Pn ~r Ps 

ip° s (-Lo + k x Us)w s = ~f a Dp' + s" P ° s Dt - gp' s + f 3 f n (U s - Un)ik x (w s - w n ) - k 2 a s z' s S(z), (16) 

Pn ' Ps Pn ~r Ps 

(k x u nt3 + k y v n>s ) = iDw nt3 , (17) 

l(-UJ + k x Un,s)p' n ,s = ~Wn,sDpl i3 , (18) 

+ k x U„, s (z = 0)]z' niS = w„, s {z = 0), (19) 

where U n ,s are the zeroth-order shear velocities (in the x direction as defined previously), u,v, and w are the perturbed 
velocities in the x, y, and z directions, p' is the first order pressure, r is the first order temperature, p is the first order 
density, and other quantities with the superscript '0' are of zeroth order. Following Chandrasekhar (1961), we eliminate all 
first-order quantities except for w n ,s- Using (11), (12), (14), and (15) we eliminate r and p' , while u n ,s and v n , s are eliminated 
using (17). We then obtain the bulk fluid equations (equations which are valid away from the interface) from (13) and (16), 
after using (18) to eliminate p n s and (19) to eliminate « n , s . 

For the bulk fluids, where all zeroth-order velocities are constant and uniform everywhere, one obtains the following 
equations which are valid throughout the superfluids and normal fluids away from the interface: 

£-D[p°{-u; + k x Us)Dw s - p° s k x (DUs)ws] - ^4-k x {U s - U n )k 2 w s 
p-' p u 

+s ) p° s D{-^ T) [-(-uj + k x Us)Dw s + k x (DU s )Ws + k x (U s - U n )DWs - k x (DU s - DUn)w s }} 
s p 

= p °(-uj + k x Un)k 2 w n - ^D[p° n (-cu + k x U n )Dw n - p° n k x {DU n )w n ] + gk 2 (Dpt) ^— - ^-k x (U s - U n )k 2 w n 

-s°p°D{-^[{-uj + k x U„)Dw n - k x (DU„)wn - k x (Us - U n )Dw n + k x (DUs - DU„)w„]} 
s p 

" fc4 -« ( ,. W l rr ) 5 M (20) 



-LJ + k x U n 



z = 



and 

p°.(-u> + k x Us)k 2 w s - ^D[p°(-u + k x U s )Dw s - p° s k x (DU s )ws] + gk 2 (Dp° s ) ^— - ^-k x (U s - U n )k 2 w s 

+ s "PsD{-^[-(-lo + k x U s )Dw s + k x (DUs)w s + k x (U s - U n )Dw s - k x (DU a - DU n )w s ]} 
s u p u 

-k A Os ( Ws . TT ) 5(z) 

\-u + k x UsJ z=0 
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= ^D[pI{-lo + k x U n )Dw n - p° n k x (DU n )w n ] - ^-k x (U s - U n )k 2 w n 

-s°p°D{-^-[{-u} + k x U n )Dw n - k x {DU n )w n - k x {U s - U n )Dw n + k x (DU s - DU n )w n ]}. (21) 
s u p u 

From (20) and (21), we see that, away from the interface, where U n , s , p n ,s, and s° are constant and uniform, w n and w s 
must satisfy 

^k x {U a - U n )(D 2 w s - k 2 w s ) = [P^Ak x (U s - U n ) - p°(-uj + k x U n )](D 2 w n - k 2 w n ), (22) 
and 

[^k x (U s - U n ) - p°(-u + k x U s )](D 2 w s - k 2 w s ) = ^k x (U s - U n )(D 2 w n - k 2 w n ). (23) 
p u p u 

One nontrivial, stratified solution is w n>s : w Ui b = N(—oj + k x U n ,B)e~ kz ,w s ,b = S(—u> + k x Us,B)e~ kz ,w„ } a = N(—cu + 
k x U n ,A)e kz , and w s ,a = S(—lu + k x lI s ,A)e kz , where the fluids on top (bottom) are labelled with subscript B (A). Another 
solution, lu = k x ll„ ± (p° s / p°) 1 ^ 2 k x (U n — U s ), describes a wave propagating in the bulk fluid; it does not describe an instability, 
so we are not concerned with it here. 

The zeroth-order quantities (except T, by assumption) jump across the interface. We integrate the equations across the 
interface to analyse the interfacial waves. Many of the terms in the above equations are of the form f(z)Dg(z), where both 
f(z) and g(z) are discontinuous functions. Direct integration of such terms is impossible in general (see Appendix A). However, 
if one assumes all the zeroth-order quantities change over the same length scale e in passing from the top fluid to the bottom 
fluid, we can integrate (20) and (21) term by term from z = — e to z = +e and then squeeze e to zero in a mathematically 
regular way. One then substitutes the solutions for (w n , s )A,B into the integrated (20) and (21) and obtains two equations 
involving the oscillation amplitudes N and S. The two equations must be compatible, i.e. linearly dependent. The condition 
for compatibility gives the following dispersion relation for waves at the interface: 

= [Ps,b(-w + k x U B>B ) 2 + p° s , A (-w + k x U a ,A) 2 + gk(p°, B ~ P° s ,a) - k 3 a s ] 
x [^(Ps, b + P°,a)(p°,b - pLa) - y(p°A + P°,B)cr n + a~f] 
+ [Pn,s(-w + k x U n ,B) 2 + p°,a(-w + k x U n ,A) 2 + gk(pZ,B - pLa) - k 3 a n ] 
x (p° n ,B + P° n ,A)(p% ^ P°,a) - y (P°,A + pIb)°s + a/3], (24) 

with 

a = ^[[( Ss + s a){p°b - P°a)(p°,b ~ P° a) + {sb - s A )(p" B + Pa)(p° : b - P°,a) 

+ (sb - saKpb - pW.,b + P° a)] + 1(sb + sa)(p% + P°a)(p% + pIa), (25) 
p = flbE-K-u; + kx U s ,b) 2 - k x {U s ,B - U n , B )(-w + k x U s , B )\ 

SBp B 

+ -^4[(-w + k x U s .A) 2 - k x {Us,A - Un,A){-UJ + k x U s , A )], (26) 
SAp A 

T = ^'1 K -1 ^ + k x U„. B ) 2 - k x (U s ,B - U n ,B)(-0J + k x U n ,B)] 

s B p B 



+ -^4[(-w + k x U n . A ) 2 - k x {U s ,A - U n , A ){-u + k x U n , A )], (27) 

SAp A 

,0 

PA,B — Ps,A,B + Pn,A,B- 

The dispersion relation (24) is a quartic in ui. It is related to the dispersion relation of the classical KH instability as 
follows. In (24), we have terms involving superfluid and normal fluid quantities, as well as new terms a/3 and 07 which involve 
the ratios of the bulk entropies. If the superfluid and normal fluid interfaces are decoupled, one must treat the interfaces 
separately using the classical analysis, i.e. excluding (3) and entropy terms in (1) and (2). For example, if either component 
is negligible, setting either p° or p° n satisfies (24) identically, or if the interfaces are allowed to oscillate separately, cross terms 
like sp s dT/dx in (1) or p s pnV na in (2) are absent. If neither component is negligible but the counterflow tends to zero, setting 
a to zero in (24) yields the classical limit, where the components essentially move together as a classical fluid and entropy 
does not affect stability. If neither component is negligible and there is significant counterflow, one must use the two-fluid 
dispersion relation (24), as there is no classical analogue. 
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2.4 Instability threshold and entropy effects 

2.4-1 Classical KH instability threshold 

In order to understand the onset of the superfluid KH instability, we first review the classical case. The KH instability occurs 
when two superimposed layers of fluids with different densities are in relative motion. When the velocity shear exceed a critical 
value, the resulting pressure gradient (from Bernoulli's Law) between the peaks and troughs of an interfacial wave overcomes 
the surface tension and gravity and the mode grows exponentially. In the classical instability, the interface oscillates at a 
frequency 



k c (p A UA + PbUb) 

PA + PB 



gk c (p A - Pb) + (k c ) 3 a (k c ) 2 pAPB{U A - U B f 



21 1/2 



pA + PB 



(pa + Pb) 2 



(28) 



for wavenumber k c (bearing in mind that the subscripts A and B denote the bottom and top fluids respectively). Given a 
particular density and shear velocity profile, the interface is unstable if Im(iii c ) > 0, which occurs for wave numbers between 



and 



fcfnax 



PaPb{Ua - Ub) 
2a(p A + Pb) 



PaPb(Ua - Ub) 2 
2a (p A + Pb) 



1 

27 



1 

+ 2~a 



P 2 aPb(Ua - Ub) 4 
(Pa + Pb) 2 



P 2 aP 2 b(Ua~U b ) 4 



(pa + Pb) 2 



- ^og(pA - ps) 



4o-g(pA - Pb) 



1/2 



1/2 



(29) 



(30) 



On the other hand, given a particular density profile, the minimum velocity shear needed to destabilise the interface depends 
on the densities of the fluids and surface tension according to (Chandrasekhar 1961) 



(U A - U B f min 



2(pA + pB) t , .,1/2 



pApB 



(31) 



2-4-2 Superfluid KH instability: a simple example 

At first glance, the superfluid dispersion relation (24) has little in common with the classical dispersion relation (28) . However, 
the underlying physics is similar, as can be seen by considering the special case p s ,A = ps,B = ps, pn,A = pn,B = pn, 
Us, a = U n ,A = 0, and U s ,b = U n ,B = U, which resembles a classical, single-fluid KH scenario. In this case, gravity cannot act 
as a stabilising force, because the densities of the top and bottom fluids are equal. However, there may easily exist another 
stabilising force which does not depend on the difference of densities (in Section 3.1.1, we discuss the possible origin of such 
a force). Thus we replace g(p s ,A — Ps,b) with F s and g(pn,A — pn,s) with F n . The dispersion relation now simplifies to 

= (2p a Lu 2 - 2p a kUio + p s k 2 U 2 - kF 3 - k 3 a 3 ) 

x - + (2 + r+ -)uj 2 - 2(1 + r)kUu + (1 + r)k 2 U 2 ] 

Pn pn r 

+ (2p„uj 2 — 2p n kUu) + p n k 2 U 2 — kF n — k [i a„) 

x [-™*± + (2 +r+-V 2 -2(1 +r)kUoj + (l + r)k 2 U 2 ], (32) 

ps ps r 

where r = sa/ sb is the entropy ratio. We can simplify this further by investigating the limit where the entropy of the bottom 
fluid is much greater than that of the top fluid, i.e. r — > oo. Our dispersion relation then becomes 

= [(2p s io 2 - 2p 3 kUcu + p 3 k 2 U 2 - kFs - k A a s ) + (2p n L0 2 - 2p n kUuj + p n k 2 U 2 - kF n - k 3 a n )](co - kU) 2 , (33) 
whose roots can be found from (u — kU) 2 = 0, which represents a stable mode, and from 

= 2(ps + p„)co 2 - 2(p s + pn)kUuj + (p s + Pn)k 2 U 2 - k(F s + F„) - k 3 (a 3 + a„), (34) 
which represents an unstable mode if 

k[2(a s + o n )k 2 - (ps + p„)U 2 k + 2(F S + F„)] < 0. (35) 
Therefore, the mode becomes unstable when U exceeds 



^min 



W(<Ts + <r n )(F s + F n ) ^ 1/4 

(ps + Pn) 2 



(36) 
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Equation (36) reduces to the classical result t/ mm = (lQFa / p 2 ) 1 ^ 4 for pA = Pb and is similar in form even for pa / Pb- 
Therefore, the underlying physics of the classical and superfluid KH instabilities are similar, except that, in the latter scenario, 
the fluids consist of coupled normal and superfluid components, only the normal fluid can transport entropy, and the top and 
bottom superfluids contain different types of Cooper pairs and hence have different entropies even at the same temperature. 
The entropy difference shifts the critical density and shear velocity as discussed below and in Figs. 1 and 2. Interestingly, we 
expect r ~ 1 under a wide range of physical conditions, yet in this regime (32) does not factorize simply into (33). To recover 
the classicahlike threshold (36), one needs r> 1, which is uncommon in nature. 

2.4-3 Superfluid KH instability: general threshold 

In the superfluid KH instability, the entropy ratio r = sa/sb shifts the instability threshold. In Fig. 1, we plot the critical 
(minimum) density ratio p = psj pA as r changes for three different values of surface tension a / ' {paUaO^ 1 ) = 0.22, 1.1, and 
17.7, while the top fluid is initially stationary. If p < p mm , the interface is stable. In Fig. 2, we plot the critical (minimum) U 
as r changes, for p = 0.1, 0.5, and 0.8. If U < t^ mm , the interface is stable. Tables 1 and 2 summarise the parameters used to 
construct these figures. From Figs. 1 and 2, we see that the shift of the instability threshold away from classical expectations 
increases as \sa — Sb\ increases. Moreover, as \sa — s B \ increases, the interface becomes unstable more easily; interfacial waves 
are induced for a wider range of k, and stability no longer depends solely on density ratio or velocity shear. Classically, the 
KH instability is driven by a difference in velocity, which causes a difference in pressure between the top and bottom layers 
of the interface by Bernoulli's law. In a superfluid, the temperature gradient Vr is also a driving force (or, more generally, 
the gradient in chemical potential). Thus, a differential shear flow at the interface drives a temperature difference between 
the layers through (3), destabilising the interface. The magnitude of this effect is related to sa/sb, since in the momentum 
equations (1) and (2) the temperature gradient force is proportional to entropy. 

In our analysis, we ignore the effects of the container wall. This is justified in the case of the free 1 So- 3 P2 interface in 
a neutron star. The friction between the fluids and the container wall also lowers the instability threshold, as concluded by 
Volovik (2002). 

3 CIRCULATION TRANSFER 

In this section, we calculate the wavelength and growth rate of the KH instability for the terrestrial Blaauwgeers et al. 
(2003) experiment and neutron star crust-core interface. In both cases, we also calculate the circulation transferred across the 
interface, which is proportional to the number of vortex lines per KH wavelength. 

The parameters pertaining to the terrestrial experiment (Blaauwgeers et al. 2003) are summarised in Table 3. The 
superfluid density is drawn from the table of total fluid densities given by Putterman (1974), assuming normal fluid densities 
between 10~ 3 p s and 10 _1 p s (the superfluid component dominated in the experiment). The surface tension at the interface 
is extrapolated from the experimental results of Osheroff & Cross (1977) and Bartkowiak et al. (2004), the latter for small 
apertures, which range from 3 x 10~ 9 to 6 x 10~ 9 N m -1 . Taking a = 10~ 8 N m _1 and applying eq. (2) from Blaauwgeers et 
al. (2003), we extrapolate the stabilising magnetic force density to be 5 N m~ 3 . 

The parameters pertaining to the neutron star crust-core interface are summarised in Table 4. Four cases are considered. 
In Case A, we investigate the instability at the crust-core interface of a typical neutron star, with density and shear velocity 
profile taken from Tsuruta (1979), Nomoto & Tsuruta (1987), Yakovlev, Levenfish & Shibanov (1999), Tsuruta (1998), and 
Svidzinsky (2003). The surface tension is extrapolated from Osheroff & Cross (1977), Bartkowiak et al. (2004), and Privorotskii 
(1975). In Cases B, C, and D, we investigate how the instability behaves as we vary the density ratio p for r = 0.1, 1, and 10. 

3.1 Blaauwgeers et al. (2003) experiment 

3.1.1 Experimental parameters 

In this laboratory experiment, He-A and He-B were superposed in differential rotation. The interface between He-A and He-B 
was magnetically stabilised and the container rotated, inducing rotation in the normal fluid on both sides of the interface, 
which then drags the superfluid component of the A-phase into rotation. The superfluid component of the B-phase was initially 
stationary. Above a critical rotational speed, the B-phase superfluid experienced sudden jumps in rotation speed, measured 
as jumps in the B-phase nuclear magnetic resonance absorption signal. This was interpreted as the transfer of vortex lines 
from the A-phase into the B-phase across the destabilized interface (Blaauwgeers et al. 2003). 

The experiment was conducted at T — 0.77T C , p = 29 bar. The A-phase normal and superfluid components as well as the 
B-phase normal fluid rotated together at 0.0039 m s _1 , while the B-phase superfluid was initially stationary. The value we 
adopt for the superfluid density, p s — 100 kg m~ 3 , is taken from the table of total fluid densities given in Putterman (1974). 
The normal fluid density is lower but an exact figure is not quoted; we try p n = 0.1, 1.0, and 10 kg m -3 (labelled as Case 1, 
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2, and 3 respectively). At the superfluid A-B interface, the order parameter of the BCS pair fluid (analogous to magnetisation 
in a fcrromagnet) changes abruptly, causing surface tension (Osheroff & Cross 1977; Bartkowiak et al. 2004). At T — 0.77T C 
and p —29 bar, experiments measure 0.3 ^ <t s /(10 _8 N m _1 ) ^ 0.6; we adopt a s — 10~ 8 N m _1 for simplicity. The magnetic 
restoring force density produced by the magnetic susceptibility gradient across the interface, F — (1/2)(xa — X-b)V(-H),) w 5 
N m -3 , where \a and \b are the susceptibilities of the A- and B-phase respectively and Kb is the magnetic field applied, is 
estimated by substituting a 3 into the approximate formula put forward by Blaauwgeers et al. (2003) (on the last page of that 
paper) to explain the observed number of vortex lines transferred across the interface. 

Gravity cannot act as the restoring force, because the top and bottom fluids are of the same densities, so the dispersion 
relation (24) needs to be modified to replace gravity with the magnetic force, viz. 

= [p°s(-w + k x U s ,B) 2 + P°a(-w + k x U a ,A) 2 - k 3 a s - kF 3 ][~(p° iA + p° B )F n - y (p°, a + P°b)oVi + 07] 

k k 3 

+ [p«,b(-w + k x U n ,B) 2 + pZ,a{~ u + k x U n ,A) 2 - k A a n - kF n ][--(p„ tA + p°, B F s - — (pZ, A + pZ,b)<t s + a/3]- (37) 

Our model treats the superfluid and normal fluid components equivalently; we are free to put the surface tension and magnetic 
restoring force into a s and F s or a n and F„. In reality, of course, there is only one surface tension and one magnetic 
restoring force, as experiments cannot distinguish between superfluid and normal fluid interfaces. However, according to the 
interpretation of the experiment tendered by Blaauwgeers et al. (2003), the A-phase and B-phase normal fluids move together 
as one continuous fluid, so we set o~ n = 0, F n = 0, a s — 10~ 8 N m -1 , and F s = 5 N m~ 3 in the absence of more detailed 
measurements. We note that the results are insensitive to this choice: setting F n = F s = 5 N m~ 3 shifts fc mm by 0.07 per 
cent for Case 1, 0.6 per cent for Case 2, and 9 per cent for Case 3. 



3.1.2 Instability threshold and growth rate 

In the Blaauwgeers et al. (2003) experiment, the normal fluid components are comoving, so the normal fluid interface is 
absolutely stable. Nevertheless, just by being there, the normal fluid significantly affects the instability threshold due to its 
ability to transport entropy via (3). This is an important result. We plot the threshold wavenumber & mm versus the entropy 
ratio r in Fig. 3. In all three cases (labelled Case 1, 2, and 3), the interface becomes unstable at fc mm different from the 
classical value (which is 7.3 x 10 3 m _1 ), particularly at lower r, where fc mm is much less than the classical value. 2 We see 
that fc mm increases with r, jumps sharply at log r ~ —0.6, then decreases gently; we find 1.0 ^ ^min/^min ^ ^ between 
—0.1 ^ logr ^ 0.1. The jump occurs at a slightly lower r in Case 3, where p n is smaller than in Cases 1 and 2. At low r, 
the fc mm is lowest (i.e. the interface is least stable) when p n is highest (Case 3), indicating that the presence of more normal 
fluid (whose interface is KH stable by assumption) helps destabilise the superfluid interface. However, this trend shifts at 
log r > 0.65, where more normal fluid density helps stabilise the interface. Interestingly, the superfluid interface is stable 
classically for, say, F s > 60 N m~ 3 but not when entropic effects are included, when it is unstable for log r < —0.6 and log 
r > 1.6. 

Classically, the normal fluid is not expected to influence the stability of the superfluid interface (in this experiment, where 
the normal fluid interface is absolutely stable, the normal fluid seems naively to act as a passive background). Yet our analysis 
shows that it does affect the stability significantly, indicating that it is driven not only by the shear flow, but also by a 'drag 
effect' (Blaauwgeers et al. 2003) between the superfluid and normal fluid components in the B-phase, the strength of which 
depends on the relative amounts of normal fluid and superfluid present. This effect comes from the term proportional to v ns 
in (1) and (2). 

In Fig. 4, we plot the growth rate of the instability [the maximum of Im(w)] for the three cases above. From Fig. 4, for 
r ~ 1, we obtain a growth time ~ 22 ms, close to the classical value. This maximum value for Im(w) is associated with cutoff 
wavenumber & cuto ff = 4.71 x 10 4 , 4.70 x 10 4 , and 4.45 x 10 4 m _1 for Cases 1, 2, and 3 respectively. Note that fc cu t ff ~ ^min- 
The instability in Case 3, where p n is highest, grows most quickly, indicating that the superfluid interface is destabilised even 
further by the normal fluid, through the drag effect. We anticipate that our prediction of the KH growth time will be tested 
by future nuclear magnetic resonance measurements. 



2 The classical value needs careful interpretation in this situation, because the Blaauwgeers ct al. (2003) experiment differs slightly from 
standard KH instability scenarios. The superfluid-supcrfluid interface resembles the classical KH setup, with the top fluid moving and 
the bottom stationary, but the normal fluid components are comoving, so the normal fluid interface is absolutely stable. Consequently, 
setting a = in (24) does not lead to the correct classical (noncntropic) limit. Instead, one should ignore the normal fluid completely 
and apply the classical formula to the superfluid-superfluid interface (Blaauwgeers et al. 2003). 
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3.1.3 Angular velocity jump 

When the interface becomes unstable, it is hypothesised that superfluid vortex lines break through the interface and transfer 
circulation from one fluid to the other. As circulation is quantized, the amount transferred depends on the wavelength of the 
KH instability, that is, on how many vortex lines (SN) can fit within one corrugation of the wave. A-phase vortex lines are 
~ fO times thicker radially than B-phase ones, so a concentration of energy is needed to push and convert A-phase lines into 
B-phase lines. This free energy is provided by the instability at the interface (Blaauwgeers et al. 2003). The precise mechanism 
of vortex transfer and conversion is still unknown, but Volovik (2002) suggested that the vortices in the A-phase are pushed 
by the Magnus force towards the vortex-free B-phase, entering the potential well formed by the corrugation of the unstable 
interface and enhancing the growth of the well until it forms a droplet of A-phase fluid filled with vorticity, which then resides 
in the B-phase (Volovik 2002). 

For Cases 1, 2, and 3, fc m j n at r « f is approximately equal to the classical value of 7300 m _1 , giving an instability 
wavelength A = 27r/fc mm = 8.6 x 10~ 4 m. The amount of circulation transferred across the interface is k = \Ua — £/b|Ao/2 
(Volovik 2002). The circulation quantum of superfluid 3 He is kq = h/2m 3 = 6.6 x 10~ 8 m 2 s" 1 , where rri3 is the mass of a 
3 He atom. Our analysis of the Blaauwgeers et al. (2003) experiment thus implies k = 1.7 x 10~ 6 m 2 s _1 and SN = k/kq w 25. 
Experimentally, Blaauwgeers et al. (2003) observed SN w f0. The discrepancy may be due, in part, to misestimating of a n ,s, 
F n , s , and/or p n<s in the absence of exact quoted values. We take heart, however, from semiquantitative agreement between 
theory and experiment and proceed to apply the theory to the crust-core interface in a neutron star. 

3.2 Neutron star 1 5 - 3 P2 interface 

3.2.1 Stellar parameters 

We assume that the phase transition from 1 So to 3 P2 neutron superfluid in a neutron star occurs at a sharp boundary near the 
crust-core interface at radius R « 6 km (Yakovlev, Levenfish & Shibanov 1999). The inner crust consists entirely of neutron 
BCS pairs in the 1 5'o state, behaving as an isotropic superfluid of density ps = 1.5 x 10 17 kg m~ 3 (Tsuruta 1979; Nomoto & 
Tsuruta 1987; Yakovlev, Levenfish & Shibanov 1999; Svidzinsky 2003). The temperature of the crust is around 10 8 K, which 
is about 0.01T C (Tsuruta 1998; Yakovlev, Levenfish & Shibanov 1999) so the normal fluid is negligible; we therefore assume 3 
p s = 1.5 x 10 17 kg m~ 3 and p n = 1.5 x 10 15 kg m 3 in the absence of detailed knowledge. The core consists of neutron BCS 
pairs in the 3 P2 state, with T = 0.2T C — 10 8 K (Yakovlev, Levenfish & Shibanov 1999). Hence the superfluid component 
is dominant here also; we assume p s = 2.8 x 10 17 kg m -3 and p n — 2.8 x 10 15 kg m -3 . This density and shear velocity 
profile we label as Case A, summarised in Table 4. We also calculate how the instability wave number fc mm changes with 
respect to density ratio p = P C rust//°core with density ratio r fixed at 0.1, 1, and 10. We present these results as Cases B, 
C, and D respectively in Table 4. Note that, in the subsonic regime we investigate, the neutron star matter is approximately 
incompressible. 

The surface tension at the neutron star crust-core interface is unknown. However, according to Privorotskii (1975), surface 
tension scales as m*ppT 2 (l — T/T c ) 3/ ' 2 , where m* is the effective mass of the particles and pp is the Fermi momentum. Noting 
that pp in a neutron superfluid is ~ 2 x 10 s times higher than terrestrial 3 He and T c is ~ 10 11 times higher, and terrestrial 
3 He is measured to have a s — 10~ 8 N m _1 (Osheroff & Cross 1977; Bartkowiak et al. 2004), we estimate that the surface 
tension of neutron star superfluid is a s ~ 10 19 N m -1 , i.e. 25 orders of magnitude higher proportionately. We are acutely 
aware of the shortcomings of this estimate and hope that more reliable calculations will be undertaken in the literature in the 
future. 

3.2.2 Glitch threshold and growth rate 

We plot ^nnn/^n-ijn versus log r for Case A in Fig. 5, where fc^ in « 0.67 m 1 is the classical value of the instability threshold, 
obtained by setting a = in (24). We see that as r increases, the interface becomes less stable. The instability wave number 
fc mm is inversely proportional to shear velocity U and directly proportional to gravitational acceleration g, as purely classical 
KH instability theory predicts (Chandrasekhar 1961). At low r, the plot is flat and fc mm is fairly close to the classical value 
(fe min = 2.662, 0.666, and 0.167 m _1 for U = 5 x 10 5 , 10 6 , and 2 x 10 6 m s _1 respectively). Case A is classically unstable, 
but the shift in fc mm is only noticeable at higher r, where the plot dips at log r > —0.5. Note that we find & mm P* 3> 1, so 
the assumption of a plane wave perturbation is justified a posteriori. 

One can also vary the density ratio p while keeping r constant. In Fig. 6, we plot fc mm against p for three different 

3 We attempted to calculate the exact temperature-dependent value of p n /p s from eqs. (3.92) and (3.94) in Vollhardt & Wolfic (2002). 
However, due to uncertainties in m* and the Fermi-liquid corrections in a neutron superfluid, we obtained p n 0.3p to ^ a j at T =0.0f T c , 
which seems to be higher than existing phenomcnological models can sustain (Link, Epstein & Baym 1993). Thus we take p n = 0.01p tota [ 
for now. 
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entropy ratios r — 0.1, 1, and 10, labelled as Cases B, C, and D respectively (see Table 4). The classical minimum wave 
number approaches infinity as the top fluid density reaches zero, i.e. the interface becomes absolutely stable as the top fluid 
disappears. When entropy is taken into account, fc mm still increases as p gets lower, but at a slower rate, especially at higher 
r. Cases B, C, and D join the classical curve at p « 0.45, 0.9, and 1.0, demonstrating again that we diverge further from the 
classical expectation as \sa — sb\ increases. 

In the foregoing, we do not consider the effects of a magnetic field and susceptibility discontinuity at the interface. In 
principle, a magnetic field can stabilise the boundary as in the Blaauwgeers experiment (see Section 3.1). We can set an 
approximate upper bound on V(flf) at B 2 /((IqR), where B is the magnetic field at the surface of the star (typically 10 s — 10 9 
T). Thus, the upper bound on the magnetic force density is between 1.06 x 10 24 and 1.06 x 10 26 N m~ 3 . However, we find 
that, for Case A (i.e. typical neutron star parameters), the minimum magnetic force density necessary to stabilise the interface 
occupies the range 2.4-2.9 x 10 38 N m -3 for —2 < logr ^ 1, much higher than the maximum force density available. 

For a given shear velocity and density profile, one can now estimate the growth rate 7 of the instability. We present 
our plot of growth rate versus r for a standard neutron star density profile (Case A) as Fig. 7; the classical growth rate is 
Im(cj c ) = 1.77 x 10 15 s _1 for these parameters. The growth rate is found to be directly proportional to the cube of the shear 
velocity and inversely proportional to surface tension. We find 7 ~ 1.017° for r ~ 1. The growth time grows monotonically 
with k and takes the value of 5.59 x 10~ 16 s at fc cu t ff = 6-53 x 10 9 m _1 (where fc cu t ff is determined by the surface tension). 

If there is no surface tension, the viscosity of the fluid suppresses the KH instability when the viscous time-scale is shorter 
that the growth time of the KH oscillation, i.e. for fc 1/ Im(w), where v is the kinematic viscosity. The viscosity of superfluid 
neutron star matter is estimated to lie between 0.1 and 0.17 m s _1 for Case A (Andersson, Comer & Glampedakis 2004). In 
Fig. 8, we present a plot of Im(w) versus r. We find Im(w) oc U 2 v~ x . For Case A, if the surface tension vanishes but viscosity 
suppresses the instability, we estimate Im(w) = l.llm(w c ) and ft cu t ff = 5.24 x 10 6 m _1 for r = 1, translating into a growth 
time of 4 x 10~ 13 s. The growth time is ~ 10 3 times larger than the limit imposed by surface tension above, indicating that 
viscosity suppresses the instability more effectively than surface tension. 

Prima facie, our KH growth time estimate seems too small; certainly, it is unresolvable by radio pulsar timing. However, 
it is highly dependent on fc cu1 - ff (^cutoff ^min> unu ke in the terrestrial experiment). As fc cu t ff depends on dissipation 
physics which is largely neglected in our model, e.g. textural terms on the right-hand sides of (1) and (2), a more complete 
analysis is needed before the KH growth time can be predicted with confidence. 

3.2.3 Angular velocity jump 

For typical neutron star parameters (Case A), near r = 1, we obtain fc m j n ~ 0.587 m _1 , and hence Ao = 10.7 m. This 
implies that the amount of circulation transferred across the interface is k — \U A - U B \\o/2 = 5.35 x 10 6 m 2 s _1 . The total 
circulation at the interface is roughly 2tyR\Ua — Ub\ ~ 6 x 10 9 m 2 s _1 . Therefore, the fractional change in angular velocity is 
Afi/fi « Ao-R ~ 1.4 x 10~ 4 . This compares surprisingly well with observations of the angular velocity jump in pulsar glitches, 
viz. 10~ 9 < AQ/Q < 10~ 4 (Shemar & Lyne 1996). 

3.2.4 Proton superfluid 

Hydrodynamic models of neutron star superfluids in the literature (Mendell 1991; Andersson, Comer & Prix 2003) often 
include protons, electrons, and muons. The latter two species can be effectively incorporated into the normal fluid in our 
analysis (Mendell 1991), but the protons cannot. How significant is the proton superfluid? 

One key effect of the protons is the drag they exert on the neutrons. The strength of the 'drag effect' between the protons 
and neutrons is set by the off-diagonal elements of the symmetrical mass density matrix, i.e. p np = p pn = Pn(m^ — m n )/m^ = 
p P (m* — m p )/m p , where p n (p p ), m n (m p ), and m* (m*) are the neutron (proton) density, bare mass, and effective mass 
respectively. Estimates of m p vary between 0.3m p and 0.7m p , while p v can be estimated from neutron star cooling rates 
(Tanigawa, Matsuzaki & Chiba 2004). If we take the proton density to be 15 per cent of the total core density of 2.8 x 10 17 kg 
m -3 , the threshold proton fraction for direct URCA cooling (Tanigawa, Matsuzaki & Chiba 2004), we obtain \p np /pn\ ~ 0.15. 
Hence the proton-neutron drag force in Eq. (14) of Mendell (1991), as a fraction of a typical term in our momentum equation 
(1), is given by (p np /p)\v (p> - v (n) ||V x v (n) |/|v (n) \u> < (p np / ' p){2kv^ ) /w) ~ 0.3 (for a typical KH instability, the phase 
velocity is of the same order as the shear velocity), where (v^) is the neutron (proton) superfluid velocity. In other 
words, the proton-neutron drag is smaller than, but comparable to, the forces we consider and should be included in general. 
It is omitted in this paper to keep the algebra manageable and to maintain focus on the entropic KH effects, presented here 
for the first time. Note that, although the normal fluid fraction under neutron star conditions may be low, the entropy ratio 
r is what matters in our analysis, and it is not negligible. 

There are other physical effects that also appear when charged species are included. In this paper, we intentionally neglect 
electromagnetic forces [see Eqs. (15) and (16) in Mendell (1991)] and electromagnetic contributions to the chemical potential 
[Eqs. (29) and (30) in Mendell (1991)], again to render the superfluid KH analysis tractable. One electromagnetic effect is 



Kelvin- Helmholtz Instability and Circulation Transfer. 



11 



included: the force term F n<s , which can stabilise the interface in principle (see Section 3.1.1 above), although in practice it 
is too weak to shift fc mm appreciably in a neutron star (see Section 3.2.2). Apart from excluding the proton superfluid, our 
hydrodynamic equations are similar to those written down by Mendell (1991). In our case, the coupling is between the normal 
fluid and superfluid components and is governed by entropy (this coupling has been unambiguously observed in laboratory 
experiments, e.g. second sound in He II). Our terms involving entropy and counterfiow are derived from the gradient of the 
chemical potential (see Section 2.1) and therefore are implicitly included in the neutron superfluid chemical potential p,„ [Eq. 
(29) in Mendell (1991)]. Another important difference is that Mendell (1991) and Andersson, Comer & Prix (2003) deal with 
a non-stratified neutron star which consists of two different but interpenetrating superfluids, whereas we deal with a stratified 
(two-zone) neutron star with two layers of superfluid in differential rotation. 



4 CONCLUSIONS 

We analyse the superfluid KH instability within the framework of two-fluid hydrodynamics including entropy flow and both 
the superfluid and normal fluid components in each fluid. We find that the coupling between the superfluid and the normal 
fluid components via entropy transport lowers the instability threshold (especially when the entropy ratio is large) even if 
there is no significant counterfiow. When the counterfiow is significant, as in the experiment conducted by Blaauwgeers et al. 
(2003), the reduction in the threshold is even greater. Classically, the KH instability occurs because the velocity shear induces 
a self-reinforcing pressure gradient across the interface. In a superfluid, the (perturbed) temperature gradient, driven by the 
differential flow, is also a driving force, the magnitude of which is controlled by ratio of the entropies of the fluids, and the 
interface becomes less stable than expected classically. 

When the interface crinkles, it is believed that superfluid vortex lines break through and transfer circulation between the 
fluids, although the mechanism is not understood at the Gross-Pitaevskii level. Since circulation is quantized, the amount 
transferred depends on how many vortex lines can fit in one corrugation of the KH surface wave. Applied to Blaauwgeers 
et al. (2003) experiment, our model yields a growth time of ~ 22 ms and threshold wavenumber fc mm = 7.3 x 10 3 m _1 , so 
that the number of vortex lines transferred is 8 N m 25. The observed number of vortex lines is 10 (Blaauwgeers et al. 2003). 
Applied to the neutron star crust-core ( 1 S - 3 P 2 ) interface, with U = 10 6 m sec and r = 1, we find Ao = 11 m, so that 
the fractional change in angular velocity is Afi/fi « XoR ~ 1.4 x 10~ 4 . This compares favourably with typical pulsar glitch 
events, which are observed to have 10~ 9 < AQ/fi < 10~ 4 (Shemar & Lyne 1996). In view of this estimate, and the high KH 
growth rate, it is possible that vortex line transfer through a surface KH instability may explain pulsar glitches. An alternative 
process, the two-stream instability of isotropic neutron and proton superfluids in differential rotation, is also an attractive 
possibility (Andersson, Comer & Prix 2003). A more thorough analysis of both, taking into account the nonuniform texture 
and anisotropic densities and thermal conductivities of the core superfluid, is worth pursuing. 
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Table 1. Fluid parameters for the calculation of threshold density ratio p = pb/pa as a function of entropy ratio r = sa/sb in 
Fig. 1. The top fluid B is stationary. The supcrfluid and normal fluid components of the bottom fluid A move at the same velocity 
V A- We investigate three values of V a (Cases I— III). The classical thresholds are P mm = 0.997, 0.947, and 0.780 for Cases I- III 
respectively. 



Parameters 


Case I 


Case II 


Case 


III 




Pn,A/pA 


9.90 x 10 -3 


9.90 x 10 -3 


9.90 


X 


10" 


■3 


Ps,a/pa 


9.90 x lO" 1 


9.90 x lO" 1 


9.90 


X 


io- 


-1 


Ub/Ua 

















Vnb-ipAUX]- 1 


17.7 


1.10 


2.18 


X 


io- 


-1 




17.7 


1.10 


2.18 


X 


io- 


-1 



Table 2. Fluid parameters for the calculation of threshold velocity shear f/ mm as a function of entropy ratio r = sa/sb in Fig. 1. 
The top fluid B is stationary. The superfluid and normal fluid components of the bottom fluid A move at the same velocity t/ mm . 
We investigate three values of p = 0.1, 0.5, and 0.8 (Cases IV- VI). 



Parameters Case IV Case V Case VI 



Pu,b/pa 


9.90 x 10" 


-4 


4.95 x 10- 


-3 


7.92 x 10- 


■3 


Ps,b/pa 


9.90 x 10- 


-2 


4.95 x 10- 


-1 


7.92 x 10- 


-1 


Pn,A/pA 


9.90 x 10" 


-3 


9.90 x 10" 


-3 


9.90 x 10" 


3 


Ps,A/pA 


9.90 x 10- 


-1 


9.90 x 10- 


-1 


9.90 x 10" 


-1 


•Tnfo-WfJiin) 4 ]- 1 


2.29 x 10- 


-3 


5.56 x 10- 


-2 


2.47 x 10- 


-1 


^b-Wt/^j 4 ]- 1 


2.29 x 10- 


-3 


5.56 x 10- 


-2 


2.47 x 10- 


-1 
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APPENDIX A: INTEGRATING OVER INTERFACIAL DISCONTINUITIES 

In general, a composite function of the form g(x)f'(x), where both f(x) and g(x) are discontinuous at an interface x = 0, 
cannot be integrated across the discontinuity. Essentially, this is because the integral J dx5(x)6(x), where S(x) and 0(x) 
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Table 3. Fluid parameters for the instability threshold calculations as a function of entropy ratio r in Figs. 3 and 4. The parameters 
are taken from the spin-up experiment of Blaauwgeers et al. (2003). The bottom and top fluids are labelled A and B respectively. 
We investigate three values of density ratio p (Cases 1-3). 



Parameters 


Units 


Case 1 


Case 2 


Case 3 


Pn,B 


kg m — 3 


0.1 


1 


10 


Ps,B 


kg m — 3 


100 


100 


100 


Pn,A 


kg m — 3 


0.1 


1 


10 


Ps,A 


kg m — 3 


100 


100 


100 


U n , B 


m s _1 


3.9 x 10~ 3 


3.9 x 10~ 3 


3.9 x 10~ 3 


U S ,B 


m s _1 


3.9 x 10~ 3 


3.9 x 10~ 3 


3.9 x 10~ 3 


U n ,A 


m s _1 


3.9 x 10~ 3 


3.9 x 10~ 3 


3.9 x 10~ 3 


U S ,A 


m s _1 











On* 


N m" 1 











<T S 


N m" 1 


1 x 10~ 8 


1 x 10~ 8 


1 x 10~ 8 




N m" 3 











F s 


N m" 3 


5 


5 


5 


uc c 
min 


m" 1 


7.3 x 10 3 


7.3 x 10 3 


7.3 x 10 3 



a Surface tensions are extrapolated from the experimental results of Osheroff & Cross (1977) and Bartkowiak et al. (2004). 
^ Magnetic restoring force is inferred from Blaauwgeers et al. (2003) and the surface tension. 
c Classical minimum instability wavenumber. 



Table 4. Fluid parameters for the instability threshold calculations as a function of entropy ratio r (Case A) and density ratio 
p (Cases B-D). The crust fluid (B) is initially stationary in all cases; the core (A). In Case A, we investigate how the threshold 
wavenumber fc mm changes with entropy ratio r = sa/sb- In Cases B, C, and D, we investigate how fc m j n changes with density 
ratio p = Pb/pa f° r r = 0.1, 1.0, and 10 while keeping the density and velocity of the core constant. 



Parameters 


Units 


Case A 




Case B 




Case C 




Case D 


9 


m s~ 2 


5 x 10 11 




5 x 10 11 




5 x 10 11 




5 x 10 11 


Pn,B 


kg m — 3 


1.5 x 10 15 


2.8 


; x 10 14 -2.8 x 10 15 


2.i 


i x 10 14 -2.8 x 10 15 


2.8 


. x 10 14 -2.8 x 10 1B 


Ps,B 


kg m -3 


1.5 x 10 17 


2.8 


; x 10 16 -2.8 x 10 17 


2.i 


i x 10 16 -2.8 x 10 17 


2.8 


; x 10 16 -2.8 x 10 17 


Pn,A 


kg m -3 


2.8 x 10 15 




2.8 x 10 15 




2.8 x 10 15 




2.8 x 10 15 


Ps,A 


kg m -3 


2.8 x 10 17 




2.8 x 10 17 




2.8 x 10 17 




2.8 x 10 17 


U n ,B 


m s _1 




















U 3 ,B 


m s _1 




















U n ,A 


m s _1 


10 6 




10 6 




10 6 




10 6 


U 3 ,A 


m s _1 


10 6 




10 6 




10 6 




10 6 


<y n a 


N m" 1 


10 19 




10 19 




10 19 




10 19 




N m- 1 


10 19 




10 19 




10 19 




10 19 


F n 


N m" 3 




















F s 


N m- 3 




















r 




0.01-10 




0.1 




1 




10 



a Surface tensions are extrapolated from Osheroff & Cross (1977), Bartkowiak ct al. (2004), and the theory in Privorotskii (1975). 

are the Dirac delta and Heaviside step functions respectively, is ill-defined. However, if one assumes that the discontinuous 
functions change smoothly over the same small length-scale e within the infinitesimal interval e < x < e, with e — > 0, the 
integral can be regularized. For example, to integrate g(x)f'(x), one writes / = ^j^x + and g = ^rx + such that 
the functions take well-defined values at the endpoints /(— e) = a, g(—e) = c, /(e) = b, and g(e) = d, yielding 

H jy X gf = ^±^. (Al) 

Another combination one encounters in (20) and (21) involves four discontinuous functions, of the form l(x)h(x)g(x)f'(x). 
One writes / = b -^x + £±2, g = ±^x + *±£, h = ^x + S±S, and I = ^x + !±l, yielding 



lim / dxlhgf 



^(b-a) [(6 + 7 )(/3 - a)(d - c) + (S - 7 )(/3 + a)(d - c) + (6 - 7 )(/? - a)(d + c)] + ±(b- a)(6 + 7 )(/3 + a)(d + c). (A2) 
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Figure 1. Instability threshold. Minimum density ratios P mm (normalised to the classical value P mm ) plotted against log r for three 
different shear velocities, as in Table 1: Case I (solid curve, lowest U, P mm = 0.997), Case II (dashed curve, intermediate U, P mm = 0.947), 
and Case III (dotted curve, greatest U, P m j n = 0.780). The interface is stable if the system lies below the curves. 




Log 10 r 

Figure 2. Instability threshold. Minimum shear velocity t/ m j n (normalised to the classical value ^ m j n plotted against log r for three 
different density ratios p, as in Table 2: Case VI (solid curve, p = 0.1), Case V (dashed curve, p = 0.5), and Case IV (dotted curve, 
p = 0.8). The interface is stable if the system lies below the curves 



One can replace the linear trial functions above with other smooth trial functions, such as trigonometric functions; the 
results do not change. We emphasize again that there is no unique result for integrals of products of discontinuous functions 
of the kind above. For example, J dx g(x)f'(x) takes different values depending on the relative length-scales over which g and 
/ change, which are determined physically. The results above pertain to the case where both g and / jump over an interval 
0(e). 



This paper has been typeset from a TfJX/ MpJX file prepared by the author. 
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Figure 3. Threshold wavenumber fc mm versus entropy ratio r for the Blaauwgeers et al. (2003) experiment. The classical threshold 
wavenumber & mm = 7.3 X 10 3 m _1 . The solid curve is Case 1 (p n ,A = Pn,B = 0.1 kg m~ 3 ), the dashed curve is Case 2 (p n ,A = Pn,B = 1 
kg m~ 3 ), and the dotted curve is Case 3 (p n ,A = Pn,B = 10 kg m~ 3 ). 




Figure 4. Maximum growth rate 7 versus entropy ratio r for the Blaauwgeers et al. (2003) experiment, normalised to the classical 
growth rate 7° = 45.10 s" 1 . The solid curve is Case 1 (p n ,A = Pn,B = 0.1 kg m~ 3 ), the dashed curve is Case 2 (p n ,A = Pn,B = 1 kg 
m~ 3 ), and the dotted curve is Case 3 {p n ,A = Pn,B = 10 kg m~ 3 ). 
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Figure 5. Threshold wavcnumber fc mm versus entropy ratio r at the crust-core boundary of a standard neutron star. This is labelled 
as Case A (t/core = 10 6 m s -1 ) is Table 4. The classical threshold wavenumber fcfL. = 0.67 m" 1 . We find that fc m ; n oc glf~ 2 . 
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Figure 6. Threshold wavenumber fc mm versus density ratio p at the crust-core boundary of a neutron star. The dashed curve is Case 
B (r = 0.1), the dashcd-dotted curve is Case C (r = 1), and the dotted curve is Case D (r = 10). 
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Figure 7. KH growth rate versus entropy ratio r for standard neutron star parameters (Case A) and zero viscosity, normalised to the 
classical growth rate (7° = 1.77 X 10 15 s _1 at U = 10 6 m s^ 1 ). The growth rate satisfies 7 oc U 3 ^^ 1 . 
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Figure 8. KH growth rate versus entropy ratio r for standard neutron star parameters (Case A), viscosity v = 0.1 m 2 s _1 , and zero 
surface tension, normalised to classical growth rate (7° = 2.27 X 10 12 s _1 ). The growth rate satisfies 7 oc U 2 v~ 1 . 



